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Abstract 



We use power-counting arguments as an organizing principle to apply chi- 
ral perturbation theory, including an explicit A, to the pp ppn^ reaction 
Q ' near threshold. There are two lowest-order leading mechanisms expected to 

^ ■ contribute to the amplitude with similar magnitudes: an impulse term, and 

a A-excitation mechanism. We examine formally sub-leading but potentially 
large mechanisms, including pion-rescattering and short-ranged contributions. 
^ ' We show that the pion-rescattering contribution is enhanced by off-shell ef- 

' fects and has a sign opposite to that of a recent estimate based on a PCAC 

pion interpolating field. Our result is that the impulse term interferes destruc- 
tively with the pion rescattering and A-excitation terms. In addition, we have 
modeled the short-ranged interaction using a and uj exchange mechanisms. 
A recoil correction to the impulse approximation is small. The total ampli- 
tude obtained including all of these processes is found to yield cross sections 
substantially smaller than the measured ones. 
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I. INTRODUCTION 



The theory of the threshold behavior of the reaction pp ppn^ has been studied for 
some time The numerical results of the early analyses were that this behavior 

is dominated by the impulse term (Fig.(|l|a)). In this process, a single pion is emitted 
from a nucleon. Effects of initial- and final-state interactions are taken into account by 
using scattering wave functions that solve the Schrodinger equation incorporating the full 
nucleon-nucleon potential. In such a process at threshold, the pion will emerge in an s-wave. 
In the early calculations the main competitor process was thought to be "pion rescattering" 
from a pion seagull term in which the nucleon emits two pions and one of them is re-absorbed 
by the other nucleon, as illustrated in Fig.(|pD). The early estimate of this pion-rescattering 
term gave a small contribution [Q] and it was believed that the impulse term alone could 
account for the essential features of the data. 

Good data for the pp — *• pp7r° reaction near threshold became available far later than 
these original theoretical estimates ■ These data clearly show that the early theoretical 
calculations were inadequate. While the impulse term can account for the energy dependence 
of the cross section (once the phase space and Coulomb effects are treated properly) the 
total cross-section is about a factor of five below the data. Clearly something is missing in 
the original theoretical analysis. To the best of our knowledge there have been two proposed 
explanations for the missing physics. 

The first explanation, due to Lee and Riska 0, is that the original analysis is lacking 
important effects due to exchanges of shorter range. In particular, they suggest that an 
effect due to a and uj meson exchanges such as that shown in Fig.(|l]c) can account for the 
discrepancy. They find that if they use meson-nucleon coupling constants and meson masses 
obtained from meson-exchange potentials this process plus the impulse approximation repro- 
duces the data within good accuracy. Their analysis was confirmed in the work of Horowitz, 
Griegel and Meyer [Q. This explanation is interesting in several ways: it is perhaps the first 
example of a shorter than pion-ranged meson-exchange effect playing such a critical role in 
an observable and it is based on a "z-graph" effect of scattering into and out of a negative 
energy state. Hence, the latter could provide further evidence for the validity of the Dirac 
phenomenology used in intermediate-energy p-nucleus collisions [0] . 

A second explanation has recently been proposed by Hernandez and Oset They 
observed that the early estimates of the pion-rescattering term were based on the on-shell 
irN scattering amplitudes, whereas in the contribution to the pp — > ppn^ reaction one of the 
pions is well off-shell. Using models based on a pion interpolating field satisfying the partially 
conserved axial current (PCAC) requirement, they estimate that the effect of the off-shell 
behavior of the amplitude in the pion-rescattering term, together with the impulse term, is 
large enough to account for the data, although they note that at present the uncertainties 
are large. 

Which of these explanations, if either, is the correct one? Clearly both effects cannot 
simultaneously be as large as the authors of the papers suggest unless there are some com- 
pensating effects. For example, if the amplitude of the impulse term were added to the 
amplitude for both of these effects, the resulting cross section would be approximately 2- 
2.5 times larger than the experimental cross section. A fundamental difficulty in trying to 
assess whether either of these explanations is correct is that almost all of the work on the 
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subject to date has been somewhat ad hoc in an important sense. The analyses have not 
been based on a systematic expansion for which one can account for all processes without 
double counting, and for which one has some a priori method of deciding which processes 
should be small. In short, there has been no simple organizing principle. 

This lack of an organizing principle has led to a reliance on intuition in deciding which 
processes are likely to be important. There are an infinite number of processes one can 
consider and the question is which processes will make significant contributions. Unfortu- 
nately, intuition can be faulty. For example, consider the process in Fig.(Pi) in which the 
interaction generates a A that emits a pion and becomes a nucleon. This process has not 
been included in any of the analyses to date. There is a simple intuitive reason why such an 
effect should be small 0: the vrA^A vertex proceeds only through p- waves in the center of 
mass of the A. At threshold, the pion is in an s-wave relative to the center of mass of the 
entire system. Thus, one can get a nonvanishing amplitude for this process only because the 
center of mass of the A may differ from the center of mass of the system, and one expects 
such differences to be small Indeed, we estimate this effect below and find it small 
compared to typical hadronic amplitudes. However, all of the effects discussed so far are 
small in this sense. When we evaluate this effect numerically its amplitude is found to be 
sensitive to the choice of the potential used to produce the scattering wave functions, being 
in some cases similar in magnitude to these other effects. 

The fact that all of these effects are small is no coincidence. Indeed, it is a straightforward 
consequence of approximate chiral symmetry. Suppose, hypothetically, that we lived in a 
chiral world in which the up- and down-quark masses were exactly zero. It is easy to establish 
the theorem that the amplitude for emission of a pion at threshold must vanish. Of course, 
in the real world the pion mass is not zero — but it also is not large compared to typical 
hadronic scales. This suggests that the amplitude can be described in a systematic manner 
as an expansion in (Q/M), where Q is a typical (small) momentum or energy scale and M 
is a typical (large- mass) hadronic scale ~ 1 GeV, such as niN, iTip or 47?/,^. This systematic 
description is called chiral perturbation theory (xPT). 

For the example just described, Q ~ characterizes the behavior of the outgoing pion 
and one has mj^/M ~ 0.1 - 0.2 . Our problem has an additional scale that is less obvious 
(although certainly well known). The internal nuclear momentum is quite large, since the 
pion mass must be generated entirely from the nuclear kinetic energy. We will show below 
that this corresponds to Q ~ (^tItt 

M)V2. This results in a larger than desirable expansion 
parameter, (m7r/M)^/^ ~ 0.4, but one that nevertheless will allow a systematic expansion. 

The simplest way to implement xPT is via effective Lagrangians aided by power-counting 
arguments. The seminal idea was contained in a paper by Weinberg |]T2|. This idea was 
developed systematically for interactions of mesons [|13| and for interactions of mesons with 
a baryon |[l^JT5[| . The generalization of these techniques to describe properties of more than 



one baryon was also due to Weinberg ||T6| and was carried out in detail in Refs. [T^JT^ • Here 
we will use Weinberg's power-counting arguments to organize a calculation of the threshold 
production of a 7r° in a collision. 

The general idea of power counting in effective Lagrangians is, in fact, far more powerful 
than xPT. In principle, one can include in the Lagrangian and the power counting any 
light degree of freedom (and not simply the Goldstone excitations). Indeed, the failure to 
treat all relevant degrees of freedom explicitly may well lead to very slow convergence for 
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the expansion and may limit its usefulness. In the case of baryons an obvious light degree 
of freedom is the A, since mA — rrij^f is small, being about 2 m^r, and the A is strongly 
coupled to the vrA^ channel. There remains some controversy over whether or not one should 
include the A explicitly. We strongly advocate the position that the A should be included 
explicitly. Moreover, we will show by explicit calculation that the leading-order amplitude 
associated with an explicit A is similar in magnitude to the impulse approximation and 
makes significant contributions. Hence, in our counting we will consider mA — to be of 
order m.,^. We note that alternative approaches assume that this quantity is of order M. For 
pion production reactions a virtual A can have a small energy denominator (~ rriTj), and it 
makes little practical sense to treat that particle as far off shell. 

The power counting has considerable power in organizing the calculation. It tell us which 
processes are expected to dominate. But it also provides strong constraints on the form of 
the low-energy Lagrangian, which in turn is reflected in constraints on effects of the off-shell 
dependence of the various sub-processes, in particular the seagull or pion rescattering. In 
principle, one ought to be able to apply the power counting in a completely systematic 
manner and work consistently to some given order. The effective theory has a finite number 
of free parameters at a given order that one fits from some set of observables. Once deter- 
mined one can predict other observables. This approach has been used with great success in 
describing the properties of the pseudoscalar mesons |13|. In the present case, for technical 
reasons which we will discuss later, this procedure cannot be implemented in a practical 
way. Instead, we will use power counting as a guide to which processes we expect to be 
large. 

This article is organized as follows: In the following section we briefly review Weinberg's 
power counting scheme in general. Next, we discuss special features of this scheme for the 
present problem. Various terms are considered. We then derive the operators associated with 
the various leading-order processes in xPT: the impulse, A-excitation, pion-rescattering, 
recoil correction, and short-range mechanism as modeled by a and u exchanges. In the 
same section these amplitudes are evaluated using three different potentials to describe the 
initial and final state pp interactions. The impulse term interferes destructively with the A- 
excitation, and pion-rescattering amplitudes. The resulting cross sections that we compute 
are found to be far smaller than the measured ones, even when the a, cu-exchange term is 
included. Next these results are discussed, and the difference between the off-shell behavior 
of our seagull term and that seen in Ref. is explained. The paper ends with a discussion 
of a possible direction for future research using xPT for pion production reactions. 

A preliminary report of these results was presented elsewhere ||19[. Before this manuscript 
was completed, a paper appeared that also uses xPT to reach the same conclusion 
regarding the interference between impulse and seagull terms. Our work differs from theirs 
in that we modify Weinberg's power counting to the particular kinematics of this problem, 
and that we consider the A explicitly. 



II. REVIEW OF WEINBERG'S POWER COUNTING 

Three-momenta Q exchanged in typical nuclear systems are on the order of the pion 
mass, Q ~ rriT^, which is small compared to the characteristic QCD mass scale: M ~ 1 
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GeV. Whenever we face such a two-scale problem it is useful to separate the corresponding 
physics by considering an effective, low-energy theory that involves only the relevant degrees 
of freedom, all with small three-momenta Q. The pion and the nucleon obviously play this 
role, since they are the lightest stable (with respect to strong interactions) hadrons. Low- 
lying resonances also ought to be relevant. Prominent among them is the A isobar of mass 
rriA = + 5, with 5 ~ 2m7r, which is comparable to Q [|14|. Moreover, consistency of 
the chiral expansion with a large- A^^c expansion requires an explicit A pi|^B|. The case 
for higher-mass baryon states is less clear, since their mass differences with respect to the 
nucleon are larger, and they couple more weakly to nucleons and pions. Higher-mass meson 
states also have masses comparable to M. In the following we will, therefore, take the view 
that the degrees of freedom that must be accounted for explicitly are the pion, the nucleon 
and the isobar, while effects of higher-mass states are included indirectly, as contributions 
to the several parameters of the effective theory. 

The low-energy parameters are not necessarily small, so the only potential expansion 
parameter is Q/M. {Q stands not only for typical three-momenta, but also for factors of 
and 6.) Note that, because we are restricted to small three-momenta, nucleons and isobars 
are non-relativistic and act very much like static sources of pions. Corrections to the static 
limit can be accounted for by an expansion in Q/mj^. 

The next task is to count powers of Q in an arbitrary diagram. This is simple to do in the 
case of diagrams with only pions [|12|, and can be straightforwardly extended to diagrams 
with one nucleon. However, systems with several nucleons require more care, due to the 
appearance of infrared enhancements in reducible graphs ||16[. Diagrams where all energy 



denominators are of order Q are called irreducible. In this case the same power counting 
applies: an irreducible diagram with Vi vertices of type i, L loops, C separately connected 
pieces, and Ef = 2 A external fermion lines, will be proportional to Q'^, with 

oo 

u = A- A + 2L-2C + J2ViAi. (1) 

i=l 

Here the so-called index of a type i vertex is defined as 

A, = rfi + I - 2, (2) 

in terms of di (which is the sum of the number of derivatives, the number of powers of m^r, 
and the number of powers of 6) and of fi (which is the number of fermion field operators). 

However, because of the non-relativistic character of nucleons, there also exist graphs 
with intermediate states that differ in energy from initial or final states by only a small 
amount of order Q'^/2m]sf. They are larger than irreducible diagrams by factors of itlm/Q- 
We call these diagrams reducible; they can be split into irreducible sub-diagrams by cutting 
only lines corresponding to initial or final particles. The sum of irreducible diagrams is just 
what is usually called the potential. Reducible diagrams can be obtained by iteration of 
the potential, generating wave functions corresponding to scattering or bound states. In 
processes such as the ones in which we are interested, the amplitude can be organized into 
an irreducible part (to which the external pions are attached) sandwiched between wave 
functions of the initial and final nuclear states, which now contain all of the reducible parts. 
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Now, if Aj > for all i, then a perturbative expansion in Q/M exists for irreducible 
diagrams. Cliiral symmetry provides exactly this constraint. All evidence suggests that 
QCD has an approximate SU{2) x SU{2) ~ S0{4) symmetry that is spontaneously broken 
to SU{2) ~ S'0(3). The pions are the Goldstone bosons associated with this breaking; as 
such, there is at least one choice of fields for which in the chiral limit (m^ 0) they couple 



via derivatives to other particles and themselves |2^, and this is sufficient to guarantee 
Ai > 0. 

The interactions of the effective low-energy Lagrangian can thus be ordered according 
to the index in Eq.(^. Below we present only those terms that are directly relevant to 
our subsequent calculation; in particular we subsume in "■ ■ ■" those interactions with addi- 
tional pion fields that are necessary to construct the correct (nonlinear) realizations of chiral 
symmetry, but are not required in the lower-order calculations that we will perform. The 



lowest-order Lagrangian is the one where Aj = for each interaction |[12|-|1^ 



= - (Vtt)^) - \mW 

+N^[ido - -^T ■ (tt X 7T)]Ar + If Art (r . ^ . V7r)Ar 

+At[2ao-5]A + |f [Art(T-5- V7r)A + /i.c.] + --- , (3) 

where = 93 MeV is the pion decay constant, 5 = — m^r is the isobar-nucleon mass 
difference, qa is the axial-vector coupling of the nucleon. Ha is the AA^vr coupling, and S 
and T are the transition spin and isospin matrices, normalized such that 

SiS^ = l{26,j - ieijkCTk) (4) 
o 

TaTf^ = -(25afe - i^abcTc)- (5) 

Notice that we defined the fields A^ and A in such a way that there is no factor of exp{—imNt) 
in their evolution. Hence rriN does not appear explicitly to this order: the baryons are static. 
We also wrote C^'^^ in the rest frame of the baryons, which is the natural choice. (Galilean 
invariance will be assured by including terms with additional derivatives.) Chiral symmetry 
determines the coefficient of the so-called Weinberg- Tomozawa term (A^^-j- . x 'k)N) but 
not of the single-pion interactions {qa, hA)- 



£«= ——[N^^^N + -^{iN^T ■ (tt X Vtt) ■ VN + h.c.)] 



The first-order Lagrangian has A, = 1 |1^|T^,[T^|T^ , 

2mjv 4/2 

+ ^N''[{c2 + cg - T^)7r^ - C3(V7r)^ - 2cim^7r^ - ^(04 + -^)eijkeabc(^krcdiiTadjTCb]N 
HNW ■7za-VN + h.c] [iN^T ■ ttS ■ V A + h.c] 



_^iVt(r ■ a ■ V7r)NN^N - ^e^^^e^,^di7TaN^a,nN N^a^r^N + ■ ■ ■ , (6) 



6 



where the q's are coefficients of 0{1/M), SrriN ~ md — rriu is the quark mass difference 
contribution to the neutron-proton mass difference, and the diS are coefficients of 0{1/ f^M). 
These seven numbers are not fixed by chiral symmetry, but it is important to point out that 
Gahlean invariance requires that the other coefficients exphcitly shown above be related to 
those appearing in C^^\ This in particular fixes the strength of the single-pion interactions 
in terms of the lowest-order coefficients qa and Ha, and of the common mass rriN. 
The second-order Lagrangian, with Aj = 2, is 

^(2) ^ d[ + ei _ _ VNN^N + h.c] 

^^HN^T ■ iraN ■ N^WN + h.c] 

2mNf^ 

+ \N^T -naxWN ■ N^aN + h.c] + ■ ■ • , (7) 

where the Cj's are other coefficients of 0{1/ f^M). 

III. POWER COUNTING AND THE PP PPn^ REACTION 

The power-counting arguments can be extended to the pp ppn^ reaction. However, 
there is one fundamental difference relative to Weinberg's standard power counting due 
to the kinematics of the present problem. In the standard power-counting arguments it 
is assumed that typical momenta carried by nucleons are ~ itlt^. However, this condition 
cannot be satisfied for the pp ppn^ reaction. Instead one finds: 

Ptyp ~ ^/m^rrrK ■ (8) 

The reason for this is quite simple. Consider the initial state in the center of mass frame. 
At threshold the total energy is 2m at + m7r, so that the initial kinetic energy of each nucleon 
is m,r/2. Since the energy is small the nonrelativistic kinetic energy formula should apply, 
p^yp/(2m7v) = w^tt/S, and Eq.(||) follows. 

While having typical momenta of order (mArm^)^/^ can alter the details of the power 
counting, it should not spoil the scheme. Indeed, the scale of the typical momenta still 
goes to zero as we approach the chiral limit and hence a chiral expansion remains sensible. 
On the other hand these momenta are characteristically larger than what one usually en- 
counters while chiral power counting and thus the expansion can be expected to be more 
slowly convergent than a more typical case. In particular the expansion parameter becomes 
y^m^/M rather than m.^/M. The point is that since a typical momentum of the nucleons 

is (mjvmTr)^^^ nothing prevents a momentum transfer of order (mArm,r)^^^ in interactions. 
Consider for example processes involving a single-pion exchange. Chiral symmetry requires 
that (to leading order), each vertex contains a derivative and hence is proportional to the 
momentum transfer. In the traditional Weinberg power counting a pion-nucleon vertex con- 
tributes (m^/M) to the total power counting of an irreducible graph. With our kinematics, 
however, those vertices would contribute {111^^711 jsfY ^ / M ~ {rriT^/MY'''^, where the second 
form follows from m^r ~ M. Similarly, in the traditional power counting a meson propa- 
gator l/(g^ — ml) contributes as while with our kinematics it goes as l/(Mm^). In 
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short, with these kinematics whenever a three momentum enters into the power counting it 
contributes to the power counting as (m^M)^/^. 

There is another subtlety in the present power-counting scheme and this concerns the 
notion of irreducibihty. It should be recalled that a sub-diagram is reducible in Weinberg's 
sense if it includes a small energy denominator ~ m^/niN and is irreducible otherwise. 
This leads to the possibility that a sub-diagram may look topologically as though it were 
reducible, in the sense that one could cut the sub-diagram into smaller sub-diagrams with 
only nucleon external legs, while in fact the diagram is irreducible. This happens if the 
kinematics requires that the energy denominator associated with cutting the sub-diagram is 
order and is not ~ m^/m^v. This is precisely what happens for the impulse term. 

There is finally a caveat. In principle, in order to implement the power counting in the 
irreducible piece, one ought to use a potential obtained from the same chiral Lagrangian 
according to the power-counting rules. While such a potential has been developed to third 
order it is only successful at energies well below the threshold for pion production. 
Accordingly, we will follow the strategy advocated by Weinberg in his treatment of the three- 
body problem [|T6| and use semi-phenomenological potentials that incorporate experimental 
information into the nucleon- nucleon interaction. This strategy has a conceptual cost: there 
may be a mismatch between the nucleon field used in the NN potential and the one based 
on the chiral Lagrangian used in our calculation of the operators. To attempt to estimate 
the scale of the uncertainty due to this we will use a number of different A^A^ potentials. 
These effects are not negligible. 

Let us now detail how this modified power counting works for the irreducible diagrams 
close to threshold, where we can restrict ourselves to s-waves. We examine each of the 
various contributions. 

Consider first the impulse term, Fig.(|^a). At first glance this appears to be order 
mj It has been recognized that the s-wave amplitude is small because in Eq.(^ 

a pion of momentum q couples to baryons only via a ■ q. Close to threshold, the interaction 
consequently proceeds via the Galilean term doTz in Eq.(|^). This yields an explicit factor 
of m^//jr, but also involves l/rriN times a gradient that contributes ptyp- Thus, the net 
contribution at the vertex is m^^"^ / (fj^M^^"^). However, this is not yet the correct order for 
this process because of the subtlety associated with the concept of irreducibihty discussed 
above. Since the outgoing pion carries an energy of the order of the pion mass, the energies 
of the A^A^ intermediate state before and after pion emission differ by ~ m^r. Therefore both 
of the intermediate states cannot simultaneously be within ~ m^/m^v of being on-shell: at 
least one intermediate state, before or after emission, is off shell by ~ m^r. 

This single, relatively-high-momentum (~ y^m^JTn/v) pion exchange must therefore be 
included in the irreducible class of operators for our process (unlike the usual case). All 
other initial- and final-state interactions will be considered reducible and included in the 
wave functions. Thus the irreducible sub-diagrams of Fig.(|l|) should in lowest order be 
drawn as in Fig.(|^). The two-nucleon interaction itself provides a factor 1//^. Indeed, 
if it originates from virtual (static) pion exchange (Fig.(0a,b)), it results from two factors 
of y^^H~m/v//7r from each vertex and one of (m^m^v)"^ from the propagator; if it arises 
from exchange of a heavier meson h (Fig.(^,d)), it is of order g%Nh/iTi1 which is typically 
~ (47r/lGeV)^ ~ ^/ fn- '^^^ inclusion of the interaction in the irreducible part also produces 
an energy denominator between the pion exchange and the pion emission. This energy 
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denominator is (-^intermediate — -E-initiai)"^ ~ kinetic energy ~ m^r. So, along with the exphcit 
factor of order w?J'^ j (/ttM^/^) from the pion emission vertex, we have an energy denominator 
of order l/m^r and a factor of 1//^ from the potential. This gives an overall contribution 
0(/-Vm./M). 

One can use similar power-counting arguments to estimate the relative contribution of 
other processes. In studying other processes we have found that they all depend on a higher 
power of the pion mass than the ^Jm^ dependence of the impulse process. Thus, power 
counting verifies the intuition that had we lived in a world with sufficiently light quarks the 
impulse term would have been dominant. Thus, the failure of the impulse term to explain 
the data by itself is interesting in that it highlights the role of chiral-symmetry breaking. 
That is, it quantifies the sense in which the up and down quark masses cannot be regarded 
as sufficiently light. 

For example, consider the A terms of Fig.(Q). The power counting here is completely 
analogous to the case just considered; the only difference is that now the intermediate 
state energy is XjiniT^ — (5), the rest-mass difference between intermediate and initial state 
being 5, while the typical kinetic energy is, as before, of order m^r. Thus, this process is 
0{f~^ \/^7r/^) or, equivalently, 0{ ^^1g^ ) relative to the impulse approximation. In 
the limit-world where the quark masses — implying 0, while — remains 

finite, A excitation would be greatly suppressed relative to the impulse term. However, in 
the real (i.e., physical) world, we have ^J"^^^ ~ 1, and this process has to be considered as 
of the same order as the impulse term. 

All other processes are suppressed compared to the leading impulse and A-excitation 
processes by powers of ^Jm-j^/M. For example, the recoil corrections to the impulse approxi- 
mation (Fig.(^)) and to the static A excitation are down by relative order mT^/M, since they 
involve two extra factors of the ratio between transferred energy (~ m^) and transferred 
momentum (~ i/rrvM). 

The pion-rescattering or seagull process of Fig.(|^) is also down by relative order m.,^/M 
compared to the leading terms. In order to see this, we first need to understand how the 
chiral power counting goes for the irirN'^N vertex. Because of its isospin structure, the 
Weinberg- Tomozawa term in Eq.(^ does not contribute. From C^^^ in Eq.(^, we see that i) 
the Ci term and ^m^r, being proportional to the cr-term, are proportional to rUg and hence 
ml; ii) the C2 + C3 term yields an interaction proportional to go-E'Tr, where go is the energy 
of the exchanged pion, so it is also oc m^; iii) the other terms, while potentially big because 
they are proportional to the momentum of the exchanged pion, contribute only to p-waves. 
The TiTiNN vertex thus contributes a factor m^/ (f^M). To complete the power counting 
for the process, note that the pion propagator is of order l/{m.j,M) and the couphng of the 
exchanged pion to the second nucleon is a/ttv^/Ztt- Combined with the seagull vertex, this 
yields a total amplitude of order 0{f~^{mT,/MY^'^). 

A few comments are in order about the seagull contribution. The first is that although 
it is formally sub-leading, it is by no means obvious that it must be very small. Dimen- 
sionless spin and isospin factors occur in most amplitudes. They can significantly enhance 
some partial waves, and these dimensionless factors are not "counted" in power-counting 
arguments. It is also clear why this contribution is likely to be much larger than estimates 
based on the on-shell vertex. The on-shell s-wave scattering amplitude is anomalously small 
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in a power-counting sense. The basic point is that on-shell the a term contribution is nearly 
equal and opposite to the C2 and C3 contributions, yielding a total much smaller than any 
of the individual terms. However, with the kinematics of pion production the terms do not 
cancel so efficiently. Thus, one expects the on-shell estimates to be low. 

Unfortunately, the recoil and the seagull processes are not the only ones of this order. 
One-loop diagrams with and without A's (e.g., those in Fig.@) also appear to this order. 
Ideally, of course, we should simply calculate these processes and include them in our anal- 
ysis. There are two difficulties with this. The first is that we are using phenomenological 
potentials and we do not how to implement the traditional xPT scheme of regulating the 
loop and absorbing the regulator dependence into the coefficient of a contact term in a way 
consistent with the potential. The second is that even if were we able to do this in a con- 
sistent way, we would not be able to make a prediction until we fit the contact term from 
some other process. Accordingly we have not included these loop diagrams in our analysis. 
We realize, of course, that these graphs could be quite significant, especially because there 
is considerable cancellation among the other terms. 

The most relevant contact terms in question are the ones in Eq.(|^), and shown in Fig.(^. 
Their finite parts are expected to be 0[^\/mT,Mm.^l {f^MrriN)) = 0(/~^(m7r/M)^/^), so it 
is also of the same order as the rescattering mechanism. In order to gauge the relevance of 
this term, we adopt here a phenomenological strategy: following Ref. [^], we assume it to be 
saturated by a and u exchanges. If this exchange model is correct, this short-range physics 
is quite significant — the amplitude is of the same magnitude as the impulse term. On the 
other hand, it is not totally obvious that the model is correct. It depends on using large 
coupling constants and a light a. At least in part, this is believed to be a phenomenological 
parameterization of some two-pion-exchange physics in the potential. It is not obvious that 
this phenomenology should be used in a context other than the potential. Accordingly we 
will first consider calculations omitting this contribution and subsequently add it in assuming 
the strength given in Ref. ||^. 

We can continue with higher-order examples, but it should be clear now how this can be 
generalized. (For simplicity we drop here the overall factor.) In the Weinberg scheme 
the amplitude associated with an irreducible diagram is of order {Q/MY where u is given 
in Eq.(|l|) and Q is the typical momentum in the problem or m,r or 6, which are taken to be 
of the same order. With the present kinematics the order of an irreducible diagram can be 
written as 

where counts the order of the explicit chiral symmetry breaking arising directly from the 
effective Lagrangian, z/g counts the order of the energy transfer in the external legs {i.e. the 
total number of time derivatives on external legs) and Vd is the dynamical order given by 
— ^ex — ^e- Thus, we cau introduce an effective index v^ff defined by 

l^eff = + Z/e + Uex = {v + Ue + Z^ex)/2 , (lO) 

and the irreducible diagram will be of order [rriTr/My^ff . 
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IV. EXPLICIT FORMS OF OPERATORS 



We now obtain the explicit forms of the various contributions by evaluating the most 
important irreducible diagrams in momentum space. Our notation is as follows: tu^ = 

+ is the energy of the (on-shell) pion produced with momentum q in the center of 
mass; p {p') is the center-of-mass momentum of the incoming (outgoing) proton labelled "1" 

(those of proton "2" are opposite); k = p — p' [hP = (p^ — p' )/2m7v) is the momentum 
(energy) transferred; cu^ = P + m^; P = p + p'; a^^^ is the spin of proton i; S = a^^'^ — a^'^^; 
and T(k) = a^^^ ■ ka^'^^ ■ k. We define the T-matrix in terms of the S-matrix via S=l+iT. 

According to the previous discussion, we expect the leading contributions to arise from 
the diagrams in Figs. @ and (|^) where the virtual pion is exchanged between static baryons. 
In the case of pion exchange with a nucleon in the intermediate state (Fig.(0a,b)), we get 



T 



IA,p 



i9A 



f.-p'T{k) -T{k)t-p 



(11) 



which is listed for comparative purposes only. We will actually calculate the impulse approx- 
imation directly from Eq.(6). Recoil corrections, as in Fig.(D, are expected to be smaller 
by a factor of m^/M, and are also included: 



T 



recoil 



tg\ {k^ 



2/3 2 2 no\2 ^-p'Tik) -T{k)E-p 

In the case of pion exchange with a delta intermediate state (Fig.(|^a,b)) 
-WaJia 1 



(12) 



18mjv/3a;i-(P)2(52-a;2 

{Pug -k-P5)t-k + ^cjg(a(^) ■ ka^^^ -Pxk- a^^^ ■ P x fca^^) . 



(13) 



Results similar to Eqs.([n]),([T^) follow for the shorter-range terms (Fig.(|^,d), Fig.(^,d)), 
and we do not write them explicitly here. The diagrams in Fig.(H) will all be included in 
Fig.(0) as far as the explicit calculation goes. Diagrams in Fig.(|^c,d) have to be accounted 
for explicitly as well as diagrams in Fig.(|^a,b). In any reasonable model, however, they 
turn out to be smaller than those in diagrams of Fig.(^,b). For example, they could arise 
from ai exchange, but then the relatively high ai mass suppresses this contribution. For the 
purpose of estimating the effect of the A, we use Eg . ([I3|) . 

There are two other corrections of order m.,^/M compared to the leading diagrams of 
Figs.(ll) and (^. Fig.(|^) represents the s-wave rescattering: 



iST 



0\2 



(C2 + C3 - 



)k°ujq — 2cim^ — 



Sm 



N 



E-k. 



(14) 



Fig.(0) is a short-range mechanism provided by Eq.(^: 



2f^mN 



(rf; + 2ei)S-P + 2ze2a 



(1) 



X a 



(2) . k 



+ 



(15) 
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Chiral symmetry tell us nothing about the strength of these terms. We can use data to 
determine the coefficients d'^, Ci, and 62- Alternatively, we can use a model to determine these 
coefficients and then try to explain the experimental results. Here we use the mechanism ffist 
proposed by Lee and Riska P] and by Horowitz et al. where the short-range interaction 
is supposed to originate from z-graphs with a and uj exchanges. In this case, 



_2^4a + CU ^(1) ^ .(2) . ^] + . . . ^ (16) 

where m„ (m^) and g„ (g^) are the mass and the (vector) coupling to nucleons of the a 
(uj) meson, and denotes the ratio of tensor to vector coupling for the u meson. (Note 
that there is a misprint in the sign of the second term of Eq. (7b) of Ref. 0; the i should 
be replaced by —i.) Comparison between Eqs.(|l^) and (p^ ) yields the model- dependent 
estimates: < + ci = -gAgl/2mlmN ^ -1.5(1/ f^M), ci = -gAgl/2mlmN ^ -2(1/ f^M), 
and 62 = gAd^i^ + C'w)/2m^mjv ~ 2(1//^M). The numerical factors are obtained from the 



Bonn OBEP A potential. Table A. 2 of reference |2^ . 

We are concerned with evaluating the matrix elements of the above operators between 
the initial ^Pq and final ^Sq pp wave functions. It is convenient to use the Reid p6|, Reid93 



27|| and Argonne V18 potentials which, for a given pp channel, are local potentials. Thus 
we evaluate the operators between coordinate space initial [i) and final (/) wave functions 
expressed by 

pr 

and 

(^1/) = ^Mr)e-'''V4^\'So) . (18) 

We convert the operators of Eqs.(12)-(|l6|) to configuration space by inverting the Fourier 
transforms. The resulting operators can then be used in configuration-space matrix elements. 
Hermiticity and [a^^^ ■ a^'^\ S ■ f] = Aia^'^^ x a^^^ ■ f lead to 



(%Ka«xa(2).f|3Po) = -(%|S-f|3Po) 

2\ 



r 



(19) 



while direct evaluation leads to 

(i5o|S-f|3Po) = -2. (20) 

We may then begin to tabulate the results. We define the matrix elements of the opera- 
tors of Eqs.(12)-(|16D as 
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M 



X 



{f\T^\z) 
PP' 



PCX) 

Jo 



(21) 



where X represents I A, recoil, etc, and H (r) is the corresponding operator, to be given 
shortly. We have extracted all of the constant factors from Eqns.(17) and (18) together with 
a sign from Eqns.(19) and (20) and placed them in front of the integral. The integral can 
be broken up into a part from zero to ro where tq is any distance greater than the range of 
the nucleon-nucleon force. For r > ro, the wave functions tio,Mi,o are linear combinations 
of regular and irregular Coulomb wave functions, so that the integral from ro to oo can be 
done analytically. Note that we can also define J-matrix elements , so that J^^ = Ji of 
Ref. [|l|. The relationship is given by 



M 



X 



iV2 



4vr gA 



PP' fn rriNm^ 



-J 



X 



(22) 



The impulse approximation is given by 



gAm^, ( d 



(23) 



The recoil term is evaluated from 



H 



recoil i 



9a 1 



' d ^ d 
4^ {Vrim^, r) - Vrim^, r)) - — 
or or 



{Vs{m^,r) - Vs{m^,r)) 



(24) 



with 



(25) 



and 



VT(/i,r) 




(26) 
(27) 



The dependence on rri^ arises from the restriction occurring at threshold that the total 
momentum of the two final nucleons be zero. In that case, we have /c° — > m.,^/2. 
The seagull term (ST) is evaluated in a similar manner from 



H 



ST I 



P 



4ci + 



6m 



N 



2ml 



- C2 + C3 - 



a: 

Imj^ 



^1 + m^r) 
47r r^ 



(28) 



One may also evaluate the graph of Fig.(^ by treating the pion-nucleon vertex in an on-shell 
approximation. In this case the term is replaced by m^^. The result is that one obtains 
an on-shell form of the seagull term (STon) with 
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H 



STon, 



9Aml 



5m M ( 

4Ci + — - - 2 C2 + C3 



9\ 



[ 2m^ V ^ ' S^^A^yj 47r 
The term involving the intermediate A is found to take the form 



9\ 



Ka \ mlm 



26 



ISiriNf^ \gA 



mt 



„MA(r) + (l-— )i?A(r) + CA(r)], 

II ijr 



(29) 



(30) 



where 

5a (r) 



3(1 + r) 



CA{r) 



An r 
2 e"™'^'' 



1 + r) ^ 2 



1 + - 



+ - 



2 1 



+ - 



m^r r 



ml 



A, 

dr ' 



(31) 



br + - • 

or r 



An r 

We also include the effects of the sigma and omega exchange terms {cr,uj). The result is 



fn rri^ 



N 



2 dr 2 dr 



where the function /^(r) accounts for exchange of the meson h between nucleons, 

AW 



Att r 



(32) 



(33) 



We follow Ref. by including the effects of form factors as defined in the Bonn potential 
23| . For completeness, we repeat that monopole form factors are used at each meson vertex 



according to the replacement 



Qh 9h 



(34) 



where A;^ is the transferred momentum and Ah is the cutoff mass. 

These expressions for H'^''^ differ from the corresponding ones of Ref. because we 
use p' = p — k instead of approximating p' by p as in the earlier work. The failure of 
this approximation has very recently been pointed out by Niskanen [^. Keeping the term 
proportional to k yields terms proportional to that are not included in Ref. However, 
this change does not effect the conclusions of Ref. , because a decrease in the contribution 
of the sigma exchange is compensated by an increase of that of the omega exchange. 

The final steps consist of computing the total matrix element Ai: 

^lA _^ ^ST _^ ^recoil _^ _yy^A ^ ^a,u> 



M 

squaring that sum, and integrating over the available phase space. We find 

rp'r^a^ dp'p'"^ q' 



(35) 



a = - 

V Jo 



rriN 



(27r)3 ' ' 2mN + uj{q')' 
where v is the laboratory velocity of the incident proton, q' is the pion momentum, uj{q') 



(36) 
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V. INPUT PARAMETERS AND RESULTS 



The various amplitudes considered in the last section depend on several parameters that 
we can determine from other processes. The impulse-approximation operator of Eq.(p3D and 



the recoil operator of Eq.(^) depend on the pion mass, m,r = 135 MeV and on 



(37) 



we use the value of Qt^nn appropriate for each potential. The A operator of Eq.(^) further 
depends on the A — N mass splitting 6 = 294 MeV [^] and on the vrA^A coupling constant, 



hA- This has been fixed from p-wave tcN scattering (see, e.g. Ref. |31 



2.1. (38) 

9A 

The seagull operator of Eq. (|28|) depends on four parameters ci,2,3 and 6mN- The q's can be 
obtained by fitting s-wave nN scattering. In Ref. ||15| they were found to be 

Cl = — 1.63/2m7v 

C2 = 6.20/2mjv (39) 
cs = — 9.86/2mAr , 

from the cr-term, the isospin-even scattering length, and the axial polarizability. Note that 
the analysis of Ref. |TB| does not include the isobar explicitly. Since the inclusion of the 
ttA^A interaction only affects s-waves at one order higher than the Cj's, the above values 
can still be used to estimate the effect of s-wave rescattering. The remaining parameter, 
Stun, can in principle also be determined from s-wave irN scattering, but would require a 
careful analysis of other isospin-violating effects. Chiral symmetry relates it to the strong 
interaction contribution to the nucleon mass splitting, which is also difficult to determine 
directly. Estimates of the electromagnetic contribution 5m jv are more reliable, 5m jy ~ — 1.5 
MeV [^, and give 5m ~ 3 MeV. To be definite, we use 

SrriN = 3 MeV. (40) 

Finally, the a, u operator of Eq.(P^ involves gu. Ah, nih, and C^, parameters listed in Table 
A. 3 of Ref. ig. 

We shall consider results of individual terms before presenting complete calculations 
using of Eq.(^). We shall often compare our results with the lUCF ^ and Uppsala |T0 



data. This allows us to understand whether or not the size of a particular term is relevant. 
The values of the cross sections divided by ■i]'^ will be displayed, where rj is the maximum 
value of the pion momentum divided by the mass of the vr". The three points at the lowest 
values of 1] are from Ref. ||10|, the remainder are from Ref. 0. We restrict our calculations 
to values of rj such that r] < 0.4, because the a ■ Vtt term of the Lagrangian (which we 
ignore) is important for higher values |^. 

We start by considering the effects of the impulse-approximation term of Eq. (|23|) . We 
compute the cross section using JH = M.^^. The results are shown in Fig.(pi), where the 
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Coulomb distortion of the proton-proton initial and final wave functions is neglected, and in 
Fig.(^), where the Coulomb effects are included. We observe that the impulse approximation 
provides cross sections much smaller than the measured ones, and that including Coulomb 
effects is important in describing the r] dependence of the new accurate cross section data. 
Henceforth, the only results we shall present include the Coulomb effects. 

The observations about the small size of the impulse term and the importance of the 
Coulomb distortion have been made before Furthermore, the small values of the 

cross sections can be traced to cancellations in the integrands. This is shown in Fig.(^0[). 
The oscillations arise from the high momentum [p ~ y/m^fTn^) of the initial state. 

Thus terms other than the impulse approximation must be included. This is also the case 
in the pp c/vr"*" reaction and in pion production on nuclear targets (e.g., see the reviews 
1 33| , |3^ ) . The computation of the terms of Eqs. (pT)) and (p^)-(|3^) is straightforward. These 
matrix elements are evaluated as a function of p' and for rj = 0.3 using the Reid potential 
(Fig.(|ni)), the Reid93 potential (Fig.(|12D) and the V18 potential (Fig.(^). The pp dTT+ 
cross sections can be understood in terms of rescattering mechanisms, so it is natural to 
include the effects of the rescattering via the seagull term. We see that including this term 
in a manner dictated by the chiral Lagrangian (Eq.(p8D) leads to a matrix element with 
a sign opposite to that of the impulse-approximation term. These two terms cancel to a 
large extent, so that one is forced to examine other terms. If one uses on-shell kinematics 
to evaluate the influence of the seagull term (Eq. ([29|) ), one finds a contribution of the same 
sign as the impulse term. 

The reason the seagull amplitude changes its sign when going from on-shell to half-off- 
shell kinematics can be seen immediately by comparing the square brackets of Eqs.(|28|) and 



(^) using the values of of Eq. (|39D and 6m ^ Eq.f^Ol). One finds for the quantity governing 
off-shell rescattering 



4ci + ^r-^ - C2 + C3 



2m£ 



2.31 

2 mjv 



(41) 



and for the quantity governing the on-shell rescattering 



4ci + -Hr - 2 C2 + C3 - 



2m1 



1.75 
2 rriN 



(42) 



Note that the contribution from isospin violation (— 5mAr/2m^) is about 10% of the total 
on-shell contribution, which is relatively large but insufficient to play a significant role in 
the current stage of understanding of the pp ppn^ reaction. 

We also discuss the other terms. The recoil term should be of order rriTt/M relative to the 
impulse term, and it is indeed quite small, as indicated in the figures. The term involving 
the intermediate A is also smaller than impulse, but depending on the potential is about 2-4 
times bigger than the recoil term. The smaller nature of the term results from cancellations 
in the integrand, as shown in Fig.(|T^. The sensitivity of the value of the matrix element to 
the choice of potential is also due to this cancellation. Qualitatively the integrands shown in 
the figure are very similar, but the net values shown in Figs.(|lT]) and (0) are very different. 

The importance of the interference between the impulse and seagull terms is underscored 
in Fig.(|l^), which shows the result of the sum of the impulse, recoil, seagull and A terms for 
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the three different potentials that we use. The correct evaluation of our Lagrangian leads 
to cross sections that are very nearly zero! Only the lowest-energy data point can be shown 
on this figure with a linear scale. Including the effects of the A in the intermediate state is 
not sufficient to provide a description of the data. None of the calculated cross sections is 
near the data, except at the lowest energy. 

This failure to account for the measured cross sections causes us to examine the effects 
of also including the effects of a and lo exchange. As shown in Fig.(|T^) we can not account 
for the cross sections with the heavy meson exchange effects. This is because the A and 
pion seagull terms interfere destructively with the other terms. 



VI. DISCUSSION 

The cross sections we compute are far smaller than the measured ones, so we cannot 
account for this process in terms of the physics we have included. However, we have learned 
something significant about the reaction. We have found processes that contribute ampli- 
tudes of a similar size to the impulse process and which should be included in any reasonable 
description: the A process, and the seagull or pion-rescattering process. The A process has 
not been considered previously for the ppyr" reaction at threshold, but it is important 

when the newest potentials are used to describe the scattering wave functions. We have 
also seen that the seagull process is highly dependent on its off-shell behavior. 

The validity of the previous conclusions is independent of chiral perturbation theory 
itself (i.e., of the convergence of the expansion in momenta). The relevance of the A isobar 
depends mostly on the observations i) that it is a relatively light degree of freedom and ii) 
that its coupling to a nucleon and an s-wave pion is essentially fixed by Galilean invariance 
in terms of its (large) coupling to a p-wave pion. The importance of the off-shell behavior 
of the s-wave vr^A^ (re) scattering results from the existence of momentum-dependent vrvr iViV 
interactions of a size similar to the momentum-independent one (which is proportional to 
the quark masses). This in turn is a consequence of chiral symmetry. 

We have also seen that these two processes interfere destructively with the impulse term. 
This is in contrast with the recent calculation of Ref. 0], where constructive interference 
was found between impulse and seagull terms. It is useful to understand the origin of this 
discrepancy. 

The key point is that the magnitude and sign of the seagull process do depend on our 
power-counting arguments. So far as the power counting is concerned, we have identified an 
interesting problem: how one proceeds when the typical momenta are ~ {m.,^mf^^^l'^ . This 
is relatively straightforward with the choices of pion and nucleon fields in which the pseudo- 
Goldstone-boson nature of the pion is manifest; in any such basis, as used above, these 
fields provide non-linear realizations of chiral symmetry. As is well-known, it is possible 
to redefine these fields without changing the physics, provided we use our field definition 
consistently throughout the calculation ||2^. Different choices of interpolating fields amount. 



after all, to no more than different bookkeeping of where various bits of the physics reside. 



For example [^, one can transform to a linear basis at the expense of producing extra 
seagull vertices that are momentum independent but not proportional to the quark masses. 
The sole purpose of such large terms is to ensure a cancellation of large pole terms and to 
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produce a small nN scattering amplitude. Such a cancellation would seem magical were it 
not understood as a consequence of a concealed chiral symmetry; power counting in this 
basis becomes highly non-trivial. Moreover, failure to include such seagull vertices in a 
calculation would spoil chiral symmetry. 

Now, the pion field in Ref. and in xPT are not equivalent. The pion interpolating field 
in Ref. is based on models exploiting PCAC whereas the pion field in the chiral Lagrangian 
is not. Thus, a priori there is nothing wrong with getting a different result for the seagull 
term contribution in the two calculations. If both calculations consistently included all other 
processes, then the differences between the present calculation and Ref. in the off-shell 
seagull contributions would be compensated for by differences in the contributions of other 
mechanisms. Indeed, starting from our general chiral Lagrangian and performing a simple 
pion-field redefinition, we can reproduce the off-shell behavior of the PCAC 'kN amplitude of 
Ref. [0; however, in doing so an extra, anomalously large contact term is generated, whose 
sole function is to cancel the effect of the change in off-shell behavior of the ttN amplitude 
brought up by the field redefinition. This is in complete analogy to the change to a linear 
basis mentioned above. Since the calculation of Ref. |0] does not include other mechanisms 
besides the impulse and seagull terms, it is rather difficult to interpret its results. 

We are not able to reproduce the data with the processes we have included, even if the 
short-ranged process is modeled by a and oj exchange using the large coupling constants 
of one boson exchange potentials. Perhaps, in face of the preceding remarks, this should 
come as no surprise. To the extent that pieces in the potential other than static OPE are 
important, we might expect similarly important contributions from one loop diagrams (such 
as the ones in Fig.(^) in the irreducible sub-diagram, again to be evaluated consistently 
with the potential. It might be no accident, then, that without such contributions we can 
not explain the data. 

The threshold behavior of pp —>■ ppn^ in principle provides us with considerable infor- 
mation about the nature of low-energy strong-interaction physics and in particular about 
the dynamics associated with approximate chiral symmetry and its spontaneous breaking. 
The amplitude for this process would be zero in the chiral limit of rrig = 0. Moreover, as 
demonstrated by the power-counting analysis in this paper, we see that for sufficiently small 
rriq the amplitude must be dominated by the impulse term. However, we have also shown 
that the chiral expansion for a process with these kinematics is essentially an expansion in 
(m7r/M)^/^. In the real world, this parameter is rather large, ~ .4, and it is by no means 
surprising that the expansion either converges slowly or fails as an asymptotic series after 
some low order. Given this large effective expansion parameter, it is not surprising that 
many processes contribute non-negligible amounts to the final amplitude. Moreover, there 
is significant cancellation between amplitudes from the various mechanisms which makes 
reliable calculations very difficult. However, even given the large uncertainties present in 
our calculations due to these cancellations, we cannot explain the data given the processes 
we have included. It is reasonably plausible that the terms we have excluded could provide a 
large enough contribution to bring our calculations into line with the observations. Clearly, 
however, developing a reliable systematic description of this processes which explains the 
data represents a large challenge. 
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FIGURES 



FIG. 1. Various contributions to the pp ppir reaction. In this and the following figures a 
single (double) solid line stands for a nucleon (Delta) and a single (double) dashed line represents 
a pion (sigma, omega); (^'/) is the wave function for the initial (final) state. 

FIG. 2. Diagrams contributing to the impulse term that are irreducible in the context of chiral 
power counting for the pp — > ppir^ reaction. 

FIG. 3. Irreducible diagrams contributing to the leading A-excitation mechanism. 

FIG. 4. Irreducible diagrams contributing to the recoil correction of the impulse term. 

FIG. 5. The seagull or pion rescattering irreducible diagram. 

FIG. 6. Some one loop irreducible diagrams. 

FIG. 7. Short-ranged irreducible diagram. 

FIG. 8. Impulse approximation. Three different potentials are used. The Coulomb distortion is 
neglected. The solid curve is obtained using the Reid potential, the dashed curve using the Reid93 
potential and the dot-dashed curve using the V18 potential. 

FIG. 9. Impulse approximation. Three different potentials are used. The Coulomb distortion is 
included. The solid curve is obtained using the Reid potential, the dashed curve using the Reid93 
potential and the dot-dashed curve using the VIS potential. 

FIG. 10. Typical integrand for the impulse term. The dashed curve is obtained using the Reid 
potential and the solid with the V18 potential. 

FIG. 11. Matrix elements as a function of p' for 77 = 0.3. See Eq.(22). The Reid potential is 
used. 

FIG. 12. Matrix elements as a function of p' for r] = 0.3. See Eq.(22). The Reid93 potential is 
used. 

FIG. 13. Matrix elements as a function of p' for rj = 0.3. See Eq.(22). The V18 potential is used. 
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FIG. 14. Integrands for the A contribution. The solid curve is obtained with the Reid and the 
dashed with the V18 potential. 

FIG. 15. M « M^"^ + M'-ecoi; ^ j^st ^ j^a 

FIG. 16. Calculations with all of our mechanisms for the Reid, Reid93, and V18 potentials, 
using eq. (32) for the heavy meson exchange effects. 
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